KDM5-mediated transcriptional activation of ribosomal protein genes alters translation efficiency to regulate mitochondrial metabolism in neurons.

Genes encoding the KDM5 family of transcriptional regulators are disrupted in individuals with intellectual disability (ID). To understand the link between KDM5 and ID, we characterized five Drosophila strains harboring missense alleles analogous to those observed in patients. These alleles disrupted neuroanatomical development, cognition and other behaviors, and displayed a transcriptional signature characterized by the downregulation of many ribosomal protein genes. A similar transcriptional profile was observed in KDM5C knockout iPSC-induced human glutamatergic neurons, suggesting an evolutionarily conserved role for KDM5 proteins in regulating this class of gene. In Drosophila, reducing KDM5 changed neuronal ribosome composition, lowered the translation efficiency of mRNAs required for mitochondrial function, and altered mitochondrial metabolism. These data highlight the cellular consequences of altered KDM5-regulated transcriptional programs that could contribute to cognitive and behavioral phenotypes. Moreover, they suggest that KDM5 may be part of a broader network of proteins that influence cognition by regulating protein synthesis.

Metacell-based differential expression analysis identifies cell type specific temporal gene response programs in COVID-19 patient PBMCs.

By profiling gene expression in individual cells, single-cell RNA-sequencing (scRNA-seq) can resolve cellular heterogeneity and cell-type gene expression dynamics. Its application to time-series samples can identify temporal gene programs active in different cell types, for example, immune cells' responses to viral infection. However, current scRNA-seq analysis has limitations. One is the low number of genes detected per cell. The second is insufficient replicates (often 1-2) due to high experimental cost. The third lies in the data analysis-treating individual cells as independent measurements leads to inflated statistics. To address these, we explore a new computational framework, specifically whether "metacells" constructed to maintain cellular heterogeneity within individual cell types (or clusters) can be used as "replicates" for increasing statistical rigor. Toward this, we applied SEACells to a time-series scRNA-seq dataset from peripheral blood mononuclear cells (PBMCs) after SARS-CoV-2 infection to construct metacells, and used them in maSigPro for quadratic regression to find significantly differentially expressed genes (DEGs) over time, followed by clustering expression velocity trends. We showed that such metacells retained greater expression variances and produced more biologically meaningful DEGs compared to either metacells generated randomly or from simple pseudobulk methods. More specifically, this approach correctly identified the known ISG15 interferon response program in almost all PBMC cell types and many DEGs enriched in the previously defined SARS-CoV-2 infection response pathway. It also uncovered additional and more cell type-specific temporal gene expression programs. Overall, our results demonstrate that the metacell-pseudoreplicate strategy could potentially overcome the limitation of 1-2 replicates.

Lineage commitment pathways epigenetically oppose oncogenic Gαq/11-YAP signaling in dormant disseminated uveal melanoma.

The mechanisms driving late relapse in uveal melanoma (UM) patients remains a medical mystery and major challenge. Clinically it is inferred that UM disseminated cancer cells (DCCs) persist asymptomatic for years-to-decades mainly in the liver before they manifest as symptomatic metastasis. Here we reveal using Gαq/11 mut /BAP wt human uveal melanoma models and human UM metastatic samples, that the neural crest lineage commitment nuclear receptor NR2F1 is a key regulator of spontaneous UM DCC dormancy in the liver. Using a quiescence reporter, RNA-seq and multiplex imaging we revealed that rare dormant UM DCCs upregulate NR2F1 expression and genes related to neural crest programs while repressing gene related to cell cycle progression. Gain and loss of function assays showed that NR2F1 silences YAP1/TEAD1 transcription downstream of Gαq/11 signaling and that NR2F1 expression can also be repressed by YAP1. YAP1 expression is repressed by NR2F1 binding to its promoter and changing the histone H3 tail activation marks to repress YAP1 transcription. In vivo CRISPR KO of NR2F1 led dormant UM DCCs to awaken and initiate relentless liver metastatic growth. Cut&Run and bulk RNA sequencing further confirmed that NR2F1 epigenetically stimulates neuron axon guidance and neural lineage programs, and it globally represses gene expression linked to G-protein signaling to drive dormancy. Pharmacological inhibition of Gαq/11 mut signaling resulted in NR2F1 upregulation and robust UM growth arrest, which was also achieved using a novel NR2F1 agonist. Our work sheds light on the molecular underpinnings of UM dormancy revealing that transcriptional programs driven by NR2F1 epigenetically short-circuit Gαq/11 signaling to its downstream target YAP1. Highlights: Quiescent solitary uveal melanoma (UM) DCCs in the liver up- and down-regulate neural crest and cell cycle progression programs, respectively.NR2F1 drives solitary UM DCC dormancy by antagonizing the Gαq/11-YAP1 pathway; small molecule Gαq/11 inhibition restores NR2F1 expression and quiescence. NR2F1 short-circuits oncogenic YAP1 and G-protein signaling via a chromatin remodeling program. Loss of function of NR2F1 in dormant UM DCCs leads to aggressive liver metastasis.

Targeted inhibition of SCFSKP2 confers anti-tumor activities resulting in a survival benefit in osteosarcoma.

Osteosarcoma(OS) is a highly aggressive bone cancer for which treatment has remained essentially unchanged for decades. Although OS is characterized by extensive genomic heterogeneity and instability, RB1 and TP53 have been shown to be the most commonly inactivated tumor suppressors in OS. We previously generated a mouse model with a double knockout (DKO) of Rb1 and Trp53 within cells of the osteoblastic lineage, which largely recapitulates human OS with nearly complete penetrance. SKP2 is a repression target of pRb and serves as a substrate recruiting subunit of the SCFSKP2 complex. In addition, SKP2 plays a central role in regulating the cell cycle by ubiquitinating and promoting the degradation of p27. We previously reported the DKOAA transgenic model, which harbored a knock-in mutation in p27 that impaired its binding to SKP2. Here, we generated a novel p53-Rb1-SKP2 triple-knockout model (TKO) to examine SKP2 function and its potential as a therapeutic target in OS. First, we observed that OS tumorigenesis was significantly delayed in TKO mice and their overall survival was markedly improved. In addition, the loss of SKP2 also promoted an apoptotic microenvironment and reduced the stemness of DKO tumors. Furthermore, we found that small-molecule inhibitors of SKP2 exhibited anti-tumor activities in vivo and in OS organoids as well as synergistic effects when combined with a standard chemotherapeutic agent. Taken together, our results suggest that SKP2 inhibitors may reduce the stemness plasticity of OS and should be leveraged as next-generation adjuvants in this cancer.

TMIGD2 is an orchestrator and therapeutic target on human acute myeloid leukemia stem cells.

Acute myeloid leukemia (AML) is initiated and sustained by a hierarchy of leukemia stem cells (LSCs), and elimination of this cell population is required for curative therapies. Here we show that transmembrane and immunoglobulin domain containing 2 (TMIGD2), a recently discovered co-stimulatory immune receptor, is aberrantly expressed by human AML cells, and can be used to identify and enrich functional LSCs. We demonstrate that TMIGD2 is required for the development and maintenance of AML and self-renewal of LSCs but is not essential for normal hematopoiesis. Mechanistically, TMIGD2 promotes proliferation, blocks myeloid differentiation and increases cell-cycle of AML cells via an ERK1/2-p90RSK-CREB signaling axis. Targeting TMIGD2 signaling with anti-TMIGD2 monoclonal antibodies attenuates LSC self-renewal and reduces leukemia burden in AML patient-derived xenograft models but has negligible effect on normal hematopoietic stem/progenitor cells. Thus, our studies reveal the function of TMIGD2 in LSCs and provide a promising therapeutic strategy for AML.

Tenascin-C expressing touch dome keratinocytes exhibit characteristics of all epidermal lineages.

The touch dome (TD) keratinocytes are specialized epidermal cells that intimately associate with the light touch sensing Merkel cells (MCs). The TD keratinocytes function as a niche for the MCs and can induce de novo hair follicles upon stimulation; however, how the TD keratinocytes are maintained during homeostasis remains unclear. scRNA-seq identified a specific TD keratinocyte marker, Tenascin-C (TNC). Lineage tracing of Tnc-expressing TD keratinocytes revealed that these cells maintain themselves as an autonomous epidermal compartment and give rise to MCs upon injury. Molecular characterization uncovered that, while the transcriptional and chromatin landscape of the TD keratinocytes is remarkably similar to that of the interfollicular epidermal keratinocytes, it also shares certain molecular signatures with the hair follicle keratinocytes. Our study highlights that the TD keratinocytes in the adult skin have molecular characteristics of keratinocytes of diverse epidermal lineages.

Functional Contribution and Clinical Implication of Cancer-Associated Fibroblasts in Glioblastoma.

PURPOSE: The abundance and biological contribution of cancer-associated fibroblasts (CAF) in glioblastoma (GBM) are poorly understood. Here, we aim to uncover its molecular signature, cellular roles, and potential tumorigenesis implications. EXPERIMENTAL DESIGN: We first applied single-cell RNA sequencing (RNA-seq) and bioinformatics analysis to identify and characterize stromal cells with CAF transcriptomic features in human GBM tumors. Then, we performed functional enrichment analysis and in vitro assays to investigate their interactions with malignant GBM cells. RESULTS: We found that CAF abundance was low but significantly correlated with tumor grade, poor clinical outcome, and activation of extracellular matrix remodeling using three large cohorts containing bulk RNA-seq data and clinical information. Proteomic analysis of a GBM-derived CAF line and its secretome revealed fibronectin (FN1) as a critical candidate factor mediating CAF functions. This was validated using in vitro cellular models, which demonstrated that CAF-conditioned media and recombinant FN1 could facilitate the migration and invasion of GBM cells. In addition, we showed that CAFs were more abundant in the mesenchymal-like state (or subtype) than in other states of GBMs. Interestingly, cell lines resembling the proneural state responded to the CAF signaling better for the migratory and invasive phenotypes. CONCLUSIONS: Overall, this study characterized the molecular features and functional impacts of CAFs in GBM, alluding to novel cell interactions mediated by CAFs in the GBM microenvironment.

SKP2 Knockout in Rb1/p53-Deficient Mouse Models of Osteosarcoma Induces Immune Infiltration and Drives a Transcriptional Program with a Favorable Prognosis.

Osteosarcoma is an aggressive bone malignancy with a poor prognosis. One putative proto-oncogene in osteosarcoma is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase. We previously demonstrated that Skp2 knockout in murine osteosarcoma improved survival and delayed tumorigenesis. Here, we performed RNA sequencing (RNA-seq) on tumors from a transgenic osteosarcoma mouse model with conditional Trp53 and Rb1 knockouts in the osteoblast lineage ("DKO": Osx1-Cre;Rb1lox/lox;p53lox/lox) and a triple-knockout model with additional Skp2 germline knockout ("TKO": Osx1-Cre;Rb1lox/lox;p53lox/lox;Skp2-/-), followed by qPCR and immunohistochemistry validation. To investigate the clinical implications of our results, we analyzed a human osteosarcoma patient cohort ("NCI-TARGET OS") with RNA-seq and clinical data. We found large differences in gene expression after SKP2 knockout. Surprisingly, we observed increased expression of genes related to immune microenvironment infiltration in TKO tumors, especially the signature genes for macrophages and to a lesser extent, T cells, B cells, and vascular cells. We also uncovered a set of relevant transcription factors that may mediate these changes. In osteosarcoma patient cohorts, high expression of genes upregulated in TKO was correlated with favorable overall survival, which was largely explained by the macrophage gene signatures. This relationship was further supported by our finding that SKP2 expression was negatively correlated with macrophage infiltration in the NCI-TARGET osteosarcoma and the TCGA Sarcoma cohorts. Overall, our findings indicate that SKP2 may mediate immune exclusion from the osteosarcoma tumor microenvironment, suggesting that SKP2 modulation in osteosarcoma may induce antitumor immune activation.

Dgcr8 functions in the secondary heart field for outflow tract and right ventricle development in mammals.

The DGCR8 gene, encoding a critical miRNA processing protein, maps within the hemizygous region in patients with 22q11.2 deletion syndrome. Most patients have malformations of the cardiac outflow tract that is derived in part from the anterior second heart field (aSHF) mesoderm. To understand the function of Dgcr8 in the aSHF, we inactivated it in mice using Mef2c-AHF-Cre. Inactivation resulted in a fully penetrant persistent truncus arteriosus and a hypoplastic right ventricle leading to lethality by E14.5. To understand the molecular mechanism for this phenotype, we performed gene expression profiling of the aSHF and the cardiac outflow tract with right ventricle in conditional null versus normal mouse littermates at stage E9.5 prior to morphology changes. We identified dysregulation of mRNA gene expression, of which some are relevant to cardiogenesis. Many pri-miRNA genes were strongly increased in expression in mutant embryos along with reduced expression of mature miRNA genes. We further examined the individual, mature miRNAs that were decreased in expression along with pri-miRNAs that were accumulated that could be direct effects due to loss of Dgcr8. Among these genes, were miR-1a, miR-133a, miR-134, miR143 and miR145a, which have known functions in heart development. These early mRNA and miRNA changes may in part, explain the first steps that lead to the resulting phenotype in Dgcr8 aSHF conditional mutant embryos.

A SOX9-B7x axis safeguards dedifferentiated tumor cells from immune surveillance to drive breast cancer progression.

How dedifferentiated stem-like tumor cells evade immunosurveillance remains poorly understood. We show that the lineage-plasticity regulator SOX9, which is upregulated in dedifferentiated tumor cells, limits the number of infiltrating T lymphocytes in premalignant lesions of mouse basal-like breast cancer. SOX9-mediated immunosuppression is required for the progression of in situ tumors to invasive carcinoma. SOX9 induces the expression of immune checkpoint B7x/B7-H4 through STAT3 activation and direct transcriptional regulation. B7x is upregulated in dedifferentiated tumor cells and protects them from immunosurveillance. B7x also protects mammary gland regeneration in immunocompetent mice. In advanced tumors, B7x targeting inhibits tumor growth and overcomes resistance to anti-PD-L1 immunotherapy. In human breast cancer, SOX9 and B7x expression are correlated and associated with reduced CD8+ T cell infiltration. This study, using mouse models, cell lines, and patient samples, identifies a dedifferentiation-associated immunosuppression mechanism and demonstrates the therapeutic potential of targeting the SOX9-B7x pathway in basal-like breast cancer.

Protective effects of macrophage-specific integrin α5 in myocardial infarction are associated with accentuated angiogenesis.

Macrophages sense changes in the extracellular matrix environment through the integrins and play a central role in regulation of the reparative response after myocardial infarction. Here we show that macrophage integrin α5 protects the infarcted heart from adverse remodeling and that the protective actions are associated with acquisition of an angiogenic macrophage phenotype. We demonstrate that myeloid cell- and macrophage-specific integrin α5 knockout mice have accentuated adverse post-infarction remodeling, accompanied by reduced angiogenesis in the infarct and border zone. Single cell RNA-sequencing identifies an angiogenic infarct macrophage population with high Itga5 expression. The angiogenic effects of integrin α5 in macrophages involve upregulation of Vascular Endothelial Growth Factor A. RNA-sequencing of the macrophage transcriptome in vivo and in vitro followed by bioinformatic analysis identifies several intracellular kinases as potential downstream targets of integrin α5. Neutralization assays demonstrate that the angiogenic actions of integrin α5-stimulated macrophages involve activation of Focal Adhesion Kinase and Phosphoinositide 3 Kinase cascades.

Molecular and network disruptions in neurodevelopment uncovered by single cell transcriptomics analysis of CHD8 heterozygous cerebral organoids.

About 100 genes have been associated with significantly increased risks of autism spectrum disorders (ASD) with an estimate of ~1000 genes that may be involved. The new challenge now is to investigate the molecular and cellular functions of these genes during neural and brain development, and then even more challenging, to link the altered molecular and cellular phenotypes to the ASD clinical manifestations. In this study, we use single cell RNA-seq analysis to study one of the top risk gene, CHD8, in cerebral organoids, which models early neural development. We identify 21 cell clusters in the organoid samples, representing non-neuronal cells, neural progenitors, and early differentiating neurons at the start of neural cell fate commitment. Comparisons of the cells with one copy of the CHD8 knockout and their isogenic controls uncover thousands of differentially expressed genes, which are enriched with function related to neural and brain development, with genes and pathways previously implicated in ASD, but surprisingly not for Schizophrenia and intellectual disability risk genes. The comparisons also find cell composition changes, indicating potential altered neural differential trajectories upon CHD8 reduction. Moreover, we find that cell-cell communications are affected in the CHD8 knockout organoids, including the interactions between neural and glial cells. Taken together, our results provide new data for understanding CHD8 functions in the early stages of neural lineage development and interaction.

Chromatin regulators in the TBX1 network confer risk for conotruncal heart defects in 22q11.2DS.

Congenital heart disease (CHD) affecting the conotruncal region of the heart, occurs in 40-50% of patients with 22q11.2 deletion syndrome (22q11.2DS). This syndrome is a rare disorder with relative genetic homogeneity that can facilitate identification of genetic modifiers. Haploinsufficiency of TBX1, encoding a T-box transcription factor, is one of the main genes responsible for the etiology of the syndrome. We suggest that genetic modifiers of conotruncal defects in patients with 22q11.2DS may be in the TBX1 gene network. To identify genetic modifiers, we analyzed rare, predicted damaging variants in whole genome sequence of 456 cases with conotruncal defects and 537 controls, with 22q11.2DS. We then performed gene set approaches and identified chromatin regulatory genes as modifiers. Chromatin genes with recurrent damaging variants include EP400, KAT6A, KMT2C, KMT2D, NSD1, CHD7 and PHF21A. In total, we identified 37 chromatin regulatory genes, that may increase risk for conotruncal heart defects in 8.5% of 22q11.2DS cases. Many of these genes were identified as risk factors for sporadic CHD in the general population. These genes are co-expressed in cardiac progenitor cells with TBX1, suggesting that they may be in the same genetic network. The genes KAT6A, KMT2C, CHD7 and EZH2, have been previously shown to genetically interact with TBX1 in mouse models. Our findings indicate that disturbance of chromatin regulatory genes impact the TBX1 gene network serving as genetic modifiers of 22q11.2DS and sporadic CHD, suggesting that there are some shared mechanisms involving the TBX1 gene network in the etiology of CHD.

Cardiac Pericytes Acquire a Fibrogenic Phenotype and Contribute to Vascular Maturation After Myocardial Infarction.

BACKGROUND: Pericytes have been implicated in tissue repair, remodeling, and fibrosis. Although the mammalian heart contains abundant pericytes, their fate and involvement in myocardial disease remains unknown. METHODS: We used NG2Dsred;PDGFRαEGFP pericyte:fibroblast dual reporter mice and inducible NG2CreER mice to study the fate and phenotypic modulation of pericytes in myocardial infarction. The transcriptomic profile of pericyte-derived cells was studied using polymerase chain reaction arrays and single-cell RNA sequencing. The role of transforming growth factor-ß (TGF-ß) signaling in regulation of pericyte phenotype was investigated in vivo using pericyte-specific TGF-ß receptor 2 knockout mice and in vitro using cultured human placental pericytes. RESULTS: In normal hearts, neuron/glial antigen 2 (NG2) and platelet-derived growth factor receptor α (PDGFRα) identified distinct nonoverlapping populations of pericytes and fibroblasts, respectively. After infarction, a population of cells expressing both pericyte and fibroblast markers emerged. Lineage tracing demonstrated that in the infarcted region, a subpopulation of pericytes exhibited transient expression of fibroblast markers. Pericyte-derived cells accounted for ~4% of PDGFRα+ infarct fibroblasts during the proliferative phase of repair. Pericyte-derived fibroblasts were overactive, expressing higher levels of extracellular matrix genes, integrins, matricellular proteins, and growth factors, when compared with fibroblasts from other cellular sources. Another subset of pericytes contributed to infarct angiogenesis by forming a mural cell coat, stabilizing infarct neovessels. Single-cell RNA sequencing showed that NG2 lineage cells diversify after infarction and exhibit increased expression of matrix genes, and a cluster with high expression of fibroblast identity markers emerges. Trajectory analysis suggested that diversification of infarct pericytes may be driven by proliferating cells. In vitro and in vivo studies identified TGF-ß as a potentially causative mediator in fibrogenic activation of infarct pericytes. However, pericyte-specific TGF-ß receptor 2 disruption had no significant effects on infarct myofibroblast infiltration and collagen deposition. Pericyte-specific TGF-ß signaling was involved in vascular maturation, mediating formation of a mural cell coat investing infarct neovessels and protecting from dilative remodeling. CONCLUSIONS: In the healing infarct, cardiac pericytes upregulate expression of fibrosis-associated genes, exhibiting matrix-synthetic and matrix-remodeling profiles. A fraction of infarct pericytes exhibits expression of fibroblast identity markers. Pericyte-specific TGF-ß signaling plays a central role in maturation of the infarct vasculature and protects from adverse dilative remodeling, but it does not modulate fibrotic remodeling.

SKP2 knockout in Rb1/p53 deficient mouse models of osteosarcoma induces immune infiltration and drives a transcriptional program with a favorable prognosis.

Purpose: Osteosarcoma (OS) is an aggressive bone malignancy with a poor prognosis. One putative proto-oncogene in OS is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase. We previously demonstrated that SKP2 knockout in murine OS improved survival and delayed tumorigenesis. Here we aim to define the SKP2 drives transcriptional program and its clinical implication in OS. Experimental Design: We performed RNA-sequencing (RNA-seq) on tumors from a transgenic OS mouse model with conditional Trp53 and Rb1 knockouts in the osteoblast lineage ("DKO": Osx1-Cre;Rb1lox/lox;p53lox/lox) and a triple-knockout model with additional Skp2 germline knockout ("TKO": Osx1-Cre;Rb1lox/lox;p53lox/lox;SKP2-/-). We validated our RNA-seq findings using qPCR and immunohistochemistry. To investigate the clinical implications of our results, we analyzed a human OS patient cohort ("NCI-TARGET OS") with RNA-seq and clinical data. Results: We found large differences in gene expression after SKP2 knockout. Strikingly, we observed increased expression of genes related to immune microenvironment infiltration in TKO tumors. We observed significant increases in signature genes for macrophages and to a lesser extent, T cells, B cells and vascular cells. We also uncovered a set of relevant transcription factors that may mediate the changes. In OS patient cohorts, high expression of genes upregulated in TKO was correlated with favorable overall survival, which was largely explained by the macrophage gene signatures. This relationship was further supported by our finding that SKP2 expression was negatively correlated with macrophage infiltration in the NCI-TARGET OS and the TCGA Sarcoma cohort. Conclusion: Our findings indicate that SKP2 may mediate immune exclusion from the OS tumor microenvironment, suggesting that SKP2 modulation in OS may induce anti-tumor immune activation.

Six3 and Six6 jointly regulate the identities and developmental trajectories of multipotent retinal progenitor cells in the mouse retina.

Formation, maintenance, and differentiation of tissue-specific progenitor cells are fundamental tasks during organogenesis. Retinal development is an excellent model for dissecting these processes; mechanisms of retinal differentiation can be harnessed for retinal regeneration toward curing blindness. Using single-cell RNA sequencing of embryonic mouse eye cups in which transcription factor Six3 was conditionally inactivated in peripheral retinas on top of germline deletion of its close paralog Six6 ("DKO"), we identified cell clusters and then inferred developmental trajectories in the integrated dataset. In control retinas, naïve retinal progenitor cells had two major trajectories leading to ciliary margin cells and retinal neurons, respectively. The ciliary margin trajectory was directly from naïve retinal progenitor cells at G1 phase, and the retinal neuron trajectory was through a neurogenic state marked by Atoh7 expression. Upon Six3 and Six6 dual deficiency, both naïve and neurogenic retinal progenitor cells were defective. Ciliary margin differentiation was enhanced, and multi-lineage retinal differentiation was disrupted. An ectopic neuronal trajectory lacking the Atoh7+ state led to ectopic neurons. Differential expression analysis not only confirmed previous phenotype studies but also identified novel candidate genes regulated by Six3/Six6 . Six3 and Six6 were jointly required for balancing the opposing gradients of the Fgf and Wnt signaling in the central-peripheral patterning of the eye cups. Taken together, we identify transcriptomes and developmental trajectories jointly regulated by Six3 and Six6, providing deeper insight into molecular mechanisms underlying early retinal differentiation.

Prerequisite endocardial-mesenchymal transition for murine cardiac trabecular angiogenesis.

Coronary heart disease damages the trabecular myocardium, and the regeneration of trabecular vessels may alleviate ischemic injury. However, the origins and developmental mechanisms of trabecular vessels remain unknown. Here, we show that murine ventricular endocardial cells generate trabecular vessels through an "angioEMT" mechanism. Time course fate mapping defined a specific wave of trabecular vascularization by ventricular endocardial cells. Single-cell transcriptomics and immunofluorescence identified a subpopulation of ventricular endocardial cells that underwent endocardial-mesenchymal transition (EMT) before these cells generated trabecular vessels. Ex vivo pharmacological activation and in vivo genetic inactivation experiments identified an EMT signal in ventricular endocardial cells involving SNAI2-TGFB2/TGFBR3, which was a prerequisite for later trabecular-vessel formation. Additional loss- and gain-of-function genetic studies showed that VEGFA-NOTCH1 signaling regulated post-EMT trabecular angiogenesis by ventricular endocardial cells. Our finding that trabecular vessels originate from ventricular endocardial cells through a two-step angioEMT mechanism could inform better regeneration medicine for coronary heart disease.

An optimized approach for multiplexing single-nuclear ATAC-seq using oligonucleotide-conjugated antibodies.

BACKGROUND: Single-cell technologies to analyze transcription and chromatin structure have been widely used in many research areas to reveal the functions and molecular properties of cells at single-cell resolution. Sample multiplexing techniques are valuable when performing single-cell analysis, reducing technical variation and permitting cost efficiencies. Several commercially available methods have been used in many scRNA-seq studies. On the other hand, while several methods have been published, multiplexing techniques for single nuclear assay for transposase-accessible chromatin (snATAC)-seq assays remain under development. We developed a simple nucleus hashing method using oligonucleotide-conjugated antibodies recognizing nuclear pore complex proteins, NuHash, to perform snATAC-seq library preparations by multiplexing. RESULTS: We performed multiplexing snATAC-seq analyses on a mixture of human and mouse cell samples (two samples, 2-plex, and four samples, 4-plex) using NuHash. The analyses on nuclei with at least 10,000 read counts showed that the demultiplexing accuracy of NuHash was high, and only ten out of 9144 nuclei (2-plex) and 150 of 12,208 nuclei (4-plex) had discordant classifications between NuHash demultiplexing and discrimination using reference genome alignments. The differential open chromatin region (OCR) analysis between female and male samples revealed that male-specific OCRs were enriched in chromosome Y (four out of nine). We also found that five female-specific OCRs (20 OCRs) were on chromosome X. A comparative analysis between snATAC-seq and deeply sequenced bulk ATAC-seq on the same samples revealed that the bulk ATAC-seq signal intensity was positively correlated with the number of cell clusters detected in snATAC-seq. Moreover, when we categorized snATAC-seq peaks based on the number of cell clusters in which the peak was present, we observed different distributions over different genomic features between the groups. This result suggests that the peak intensities of bulk ATAC-seq can be used to identify different types of functional loci. CONCLUSIONS: Our multiplexing method using oligo-conjugated anti-nuclear pore complex proteins, NuHash, permits high-accuracy demultiplexing of samples. The NuHash protocol is straightforward, works on frozen samples, and requires no modifications for snATAC-seq library preparation.

Regulation of ERα-dependent breast cancer metastasis by a miR-29a signaling.

Malignant breast cancer (BC) remains incurable mainly due to the cancer cell metastasis, which is mostly related to the status of Estrogen receptor alpha (ERα). However, our understanding of the mechanisms through which ERα regulates cancer cell metastasis remains limited. Here we identified a miR-29a-PTEN-AKT axis as a downstream signaling pathway of ERα governing breast cancer progression and metastasis. Two estrogen response element (ERE) half sites were identified in the promoter and enhancer regions of miR-29a, which mediated transcriptional regulation of miR-29a by ERα. Low level of miR-29a showed association with reduced metastasis and better survival in ERα+ luminal subtype of BC. In contrast, high level of miR-29a was detected in ERα- triple negative breast cancer (TNBC) in association with distant metastasis and poor survival. miR-29a overexpression in BC tumors increased the number of circulating tumor cells and promoted lung metastasis in mice. Targeted knockdown of miR-29a in TNBC cells in vitro or administration of a nanotechnology-based anti-miR-29a delivery in TNBC tumor-bearing mice in vivo suppressed cellular invasion, EMT and lung metastasis. PTEN was identified as a direct target of miR-29a, inducing EMT and metastasis via AKT signaling. A small molecular inhibitor of AKT attenuated miR-29a-induced EMT. These findings demonstrate a novel mechanism responsible for ERα-regulated breast cancer metastasis, and reveal the combination of ERα status and miR-29a levels as a new risk indicator in BC.